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Over  fifty  years  ago,  Carl  Eckart  described  the  sequence 
of  processes  as  a passive  scalar  is  stirred  and  mixed  in  a 
turbulent  flow  {Eckart,  1948).  At  first,  during  the  stirring 
phase,  the  variance  of  the  scalar  gradient  is  greatly  in- 
creased, but  later,  during  the  mixing  phase,  the  gradients 
become  sufficiently  sharp  that  molecular  diffusion  be- 
comes important  and  the  gradient  variance  rapidly  de- 
creases. The  process  is  of  great  practical  importance  in 
many  engineering  situations  as  well  as  being  familiar  from 
adding  cream  to  coffee. 

The  interplay  of  stirring  and  mixing  is  also  important  in 
the  ocean  for  temperature,  salinity,  chlorophyll  and  other 
naturally  occurring  tracers  as  well  as  for  introduced  mate- 
rial. The  situation  is  complicated  by  the  numerous  physical 
processes  that  can  cause  stirring  and  by  the  interplay  of 
processes  that  act  along  isopycnals  (surfaces  of  constant 
mean  density,  where  “mean”  requires  careful  definition) 
with  “diapycnal”  processes  acting  across  the  density  sur- 
faces. 

The  Twelfth  ‘Aha  Huliko’a  Hawaiian  Winter  Workshop 
was  convened  from  15  to  19  January  2001  to  review  this 
theme,  with  the  support  of  the  Office  of  Naval  Research 
and  the  participation  of  25  invited  speakers.  The  general 
topics  included  in  presentations  and  discussions  were  (1) 
mesoscale  eddy  parameterization,  (2)  general  approaches 
to  stirring,  including  the  application  of  ideas  from  dynami- 
cal systems  theory,  (3)  inertial  instability,  submesoscale 
motions  and  vortical  motions,  (4)  the  interplay  of  isopyc- 
nal  and  diapycnal  processes,  (5)  processes  in  the  surface 
mixed  layer,  (6)  the  stirring  and  mixing  of  active,  particu- 
larly biological,  tracers,  (7)  mixing  efficiency,  i.e.  the  frac- 
tion of  energy  lost  from  the  mean  flow  which  produces  a 
vertical  buoyancy  flux,  and  (8)  differential  mixing  of  heat 
and  salt  in  doubly-stable  environments  as  well  as  those 
which  permit  double  diffusion.  Before  summarizing  the 
discussions  on  these  topics,  we  first  review  the  general 
approach  to  parameterizing  unresolved  processes. 

The  need  for  parameterization 

Ocean  general  circulation  models  used  in  climate  stud- 
ies lack  the  resolution  to  treat  small-scale  processes  explic- 
itly and  so  must  parameterize  their  effects  in  terms  of  the 


average  properties  that  are  resolved.  If  the  eddy  flux  of  a 
scalar  is  related  by  a tensor  to  the  local  gradient  of  the 
mean  concentration  of  the  scalar,  as  for  a mixing  length 
argument,  the  symmetric  part  of  the  tensor  is  diagonaliz- 
able  and  may  be  thought  of  as  anisotropic  diffusion.  The 
general  assumption  is  that  two  principal  directions  are 
along  the  “neutral”  surfaces  of  constant  mean  potential 
density,  with  very  much  greater  diffusivity  in  these  direc- 
tions than  in  the  direction  normal  to  the  surfaces.  The 
isopycnal  diffusion  along  the  neutral  surfaces  is  largely 
associated  with  mesoscale  eddies,  whereas  the  diapycnal 
mixing  across  them  is  caused  by  breaking  internal  waves 
and,  in  some  locations,  double-diffusive  processes. 

The  antisymmetric  part  of  the  tensor  describes  a “skew 
flux”  normal  to  the  mean  gradient  of  the  scalar.  Part  of  this 
flux  is  non-divergent  and  does  not  affect  the  mean  scalar 
concentration.  The  remaining,  divergent,  part  may  be  ex- 
pressed as  simple  advection  of  the  scalar  by  an  additional 
mean  flow,  the  “skew  velocity”,  over  and  above  that  which 
would  be  recorded  by  current  meters  at  a fixed  location.  To 
some  extent  the  skew  flux  can  be  thought  of  as  represent- 
ing an  eddy  flux  of  thickness,  or  spacing  between 
neighbouring  isopycnals,  but  it  is  actually  a more  subtle 
representation,  in  a coordinate  system  using  z as  the  verti- 
cal coordinate,  of  a flux  that  would  appear  in  a coordinate 
system  using  the  potential  density  instead. 

These  results  are  purely  kinematic.  Specific  numerical 
values  for  the  diffusivities  and  skew  velocities  could  be 
obtained  from  suitable  measurements  in  the  present  ocean. 
Any  model  that  aspires  to  predictive  capability,  however, 
requires  an  understanding  of  the  dynamical  processes  re- 
sponsible for  the  transports  to  which  the  model  output  is 
sensitive.  The  parameterizations  must  then  be  specified  as 
formulae  rather  than  just  numerical  values.  Consider 
diapycnal  mixing  by  internal  wave  breaking.  If  it  is  to  be 
characterized  not  by  a number  but  by  a formula  then  the 
origin  of  the  wave  field  becomes  important.  Vastly  differ- 
ent feedback  loops  can  be  obtained  in  climate  studies  de- 
pending on  whether  the  intensity  of  the  wave  field  and 
hence  diapycnal  mixing  is  determined  by  the  ocean  circu- 
lation, cascading  energy  via  eddies  to  the  waves,  or  by  the 
atmospheric  wind  field  generating  waves  at  the  sea  sur- 
face, or  by  the  surface  tide  generating  waves  at  the  sea 
floor. 
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The  need  to  parameterize  subgrid  processes  is  not  idio- 
syncratic to  oceanic  general  circulation  models  but  generic 
to  all  models.  There  are  always  unresolved  processes. 

Mesoscale  eddies 

So-called  mesoscale  eddies,  with  a spatial  scale  of  the 
order  of  1 00  km  and  a time  scale  of  the  order  of  a month, 
are  a dominant  feature  of  ocean  variability.  They  are  dy- 
namically analogous  to  synoptic-scale  eddies  in  the  atmos- 
phere. The  deformation  of  an  air  layer  by  such  eddies  is 
depicted  in  Pierre  Welander’s  (1955)  now  famous  check- 
erboard figure,  here  reproduced  as  Figure  1.  The  figure 
also  depicts  the  stirring  along  oceanic  isopycnals  if  time 
and  space  are  appropriately  rescaled. 
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Figure  1.  Stirring  without  mixing  of  passive  tracer.  The  initial 
streamlines  of  the  two-dimensional  turbulent  flow  are  shown  on 
top.  The  flow  deforms  the  checkerboard  pattern  into  elongated 
filaments.  In  the  atmospheric  case,  the  length  of  the  checkerboard 
is  1200  km  and  the  deformation  is  shown  after  6,  12,  24  and  36  h. 
In  the  oceanic  case,  the  length  of  the  checkerboard  is  1 00  km  and 
the  deformation  is  shown  after  3,  6,  12  and  18  d.  (After  We- 
lander,  1955). 

Determining  the  magnitudes  of  the  isopycnal  diffusivity 
and  skew  velocity  associated  with  oceanic  mesoscale  ed- 
dies requires  an  understanding  derived  from  a combination 
of  general  principles,  dynamical  studies  of  processes  such 
as  baroclinic  instability,  and  numerical  models  that  resolve 
the  eddies  in  various  idealized  situations.  There  are  no 
simple  answers.  Invoking  general  principles,  such  as  mix- 
ing potential  vorticity,  can  imply  unphysical  sources  of 


momentum  and  energy;  eddy  parameterizations  should  not 
do  what  eddies  cannot  do!  To  some  extent  the  eddies 
serve  to  release  available  potential  energy  from  the  mean 
state,  as  in  the  popular  scheme  of  Gent  and  McWilliams 
(1990),  but  this  will  lead  to  an  ocean  at  rest  in  a situation, 
without  forcing  or  dissipation,  where  a steady,  stable,  cir- 
culation is  more  likely.  Moreover,  idealized  eddy- 
resolving runs  have  mostly  been  used  to  produce  numeri- 
cal values  for  the  free  parameter  in  the  Gent  and 
McWilliams  scheme,  rather  than  the  formula  in  terms  of 
resolved  variables  which  we  really  need  for  predictive 
models. 

The  extant  parameterizations  of  mesoscale  eddies  really 
just  treat  the  stirring  effect  of  the  eddies.  One  unanswered 
question  is  whether  this  stirring  connects  directly  to  mo- 
lecular mixing  at  very  fine  scales,  or  whether  there  is  an 
intermediate  mechanism,  in  series  with  the  eddy  stirring, 
which  accomplishes  the  transfer  to  mixing.  Another  ques- 
tion concerns  the  ultimate  fate  of  the  energy  released  from 
the  mean  state.  Parameterization  schemes  essentially  as- 
sume that  this  is  dissipated  adiabatically,  but  this  may  not 
be  true.  We  return  to  both  these  issues  later. 

Another  very  serious  question  is  whether  the  effect  of 
mesoscale  eddies  can  be  represented  as  a purely  local  ef- 
fect, with  eddy  fluxes  just  proportional  to  local  mean  gra- 
dients. This  may  be  reasonable  in  some  regions,  such  as 
the  Antarctic  Circumpolar  Current,  but  inappropriate  for 
non-zonal  flows  such  as  the  Gulf  Stream  from  which  radi- 
ating Rossby  waves  and  coherent  rings  may  carry  proper- 
ties over  large  enough  distances  to  invalidate  any  local 
mixing  length  approach.  In  such  situations  eddies  may 
drive,  rather  than  dissipate,  the  mean  flow. 

Ocean  models  may  remain  inadequate  until  mesoscale 
eddies  are  resolved  explicitly.  Quite  apart  from  the  diffi- 
culties associated  with  mesoscale  eddies,  high  resolution 
also  seems  to  be  required  to  produce  a sufficiently  high 
effective  Reynolds  number  for  the  correct  representation  of 
phenomena  such  as  the  separation,  variability,  and  penetra- 
tion of  the  Gulf  Stream.  Although  results  stabilize  with 
increasing  resolution,  they  remain  sensitive  to  choices  in 
forcing  and  viscosity  parameterization,  even  at  the  finest 
grid  spacing  used  today. 

General  approaches  to  stirring 

Turbulent  stirring  leads  to  the  distortion  of  a tracer  field. 
The  exact  distortion  is  determined  by  the  complete  La- 
grangian  evolution  of  the  flow  field,  but  much  can  be  in- 
ferred from  the  velocity  gradient  tensor.  If  the  relative  mo- 
tion near  a point  is  dominated  by  vorticity,  then  tracer  gra- 
dients tend  to  align  with  the  shear  and  not  increase.  If,  on 
the  other  hand,  the  rate  of  strain  dominates  the  relative 
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motion,  then  a tracer  tends  to  become  elongated  along  the 
principal  axis  having  the  largest  rate  of  strain  and  the 
tracer  gradient  grows  exponentially. 

Simulations  of  two-dimensional  flows  show  additional 
phenomena  such  as  barriers  to  mixing  and  areas  of  fila- 
ment production.  These  barriers  and  areas  can  be  identified 
when  additional  information  is  used.  In  one  approach,  the 
vorticity  (from  the  velocity  field)  is  combined  with  the 
rotation  of  the  rate  of  strain  tensor  (from  the  acceleration 
field)  into  an  effective  rotation  to  characterize  the  flow 
field.  Another  approach  uses  concepts  from  dynamical 
systems  theory.  It  identifies  hyperbolic  regions  in  the  flow 
field  (characterized  by  strong  deformations  in  the  velocity 
gradient  tensor)  to  calculate  special  material  curves,  or 
manifolds,  that  determine  the  stirring  and  mixing.  A blob 
of  tracer  can  advect  without  much  distortion,  but  then  be 
ripped  apart  rather  abruptly  when  it  encounters  these  hy- 
perbolic regions  in  a flow.  Realistic  simulations  of  the  sur- 
face currents  in  the  Gulf  of  Mexico  illustrate  this  process 
(Figure  2). 


Figure  2.  Evolution  of  a blob  of  tracer  in  a simulation  of  the 
surface  currents  in  the  Gulf  of  Mexico.  The  Loop  Current  is  lo- 
cated in  the  lower  right  comer  and  a cyclone  to  the  northwest. 
Height  anomalies  are  shown  as  dashed  lines.  The  red  and  blue 
lines  are  inflowing  and  outflowing  manifolds.  There  are  two  hy- 
perbolic regions  associated  with  the  cyclone  and  the  Loop  Cur- 
rent. On  day  196.75  a circular  blob  (black)  straddles  an  inflowing 
manifold  (blue).  This  manifold  delineates  which  portion  of  the 
blob  will  travel  west  into  the  central  Gulf  and  which  portion  will 
travel  east  through  the  Straits  of  Florida.  The  blob  is  stirred  with 
modest  deformation  for  the  next  3.5  days  as  it  rotates  around  the 
cyclone.  On  day  200.25,  however,  the  blob  approaches  the  hy- 
perbolic point  created  by  the  intersecting  inflowing  and  outflow- 
ing (blue)  manifolds  to  the  southwest  of  the  cyclone.  Here  it  be- 
gins to  show  strong  compression  along  the  inflowing  manifold 
direction  and  strong  distortion  along  the  outflowing  manifold 
direction.  Day  202.00  shows  the  fate  of  the  blob  when  it  reaches 
this  intersection.  It  has  collapsed  along  the  outflowing  manifold. 
(Courtesy  of  A.D.  Kirwan  and  L.  Kantha). 
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Inertial  instability,  submesoscale  motions,  and 
vortical  motions 

Mesoscale  motions  may  arise  from  baroclinic  instability 
of  large-scale  flows  in  the  ocean  and  may,  in  turn,  be  un- 
stable to  submesoscale  eddies  an  order  of  magnitude 
smaller.  Near  the  sea  surface  it  seems  that  these  eddies, 
often  visible  at  the  sea  surface  via  beautiful  spiral  patterns 
(Figure  3),  are  mainly  cyclonic,  quite  possibly  because 
anticyclonic  ones  would  break  down  rapidly  via  an  inertial, 
or  barotropic,  instability.  These  small  eddies  may  be  a link 
in  the  pathway  from  stirring  to  mixing,  but  this  has  not 
been  quantified.  Similarly,  in  the  ocean  interior  there 
seems  to  be  a family  of  “vortical  modes”,  of  largely  un- 
known scale  range  and  intensity,  which  play  a role  in  lat- 
eral stirring.  They  can  be  generated  by  frictional  effects  in 
flow  past  topographic  features  and,  possibly,  by  the  col- 
lapse and  geostrophic  adjustment  of  regions  mixed  by 
breaking  internal  waves. 


Figure  3.  A pair  of  interconnected  spirals  in  the  Mediterranean 
Sea  south  of  Crete.  The  vortex  pair  has  a clearly  visible  stagna- 
tion point  between  the  two  spirals,  the  cores  of  which  are  aligned 
with  the  preconditioning  wind  field.  Picture  taken  on  October  7, 
1984.  (Courtesy  of  W.  Munk). 

A general  question  about  the  links  between  low  fre- 
quency eddies  and  higher  frequency  motions  concerns  the 
extent  to  which  the  mesoscale  or  submesoscale  eddies 
generate  internal  waves,  particularly  in  situations  where 
the  eddies  generate  sharp  fronts.  Theoretically,  the 
mesoscale  and  submesoscale  eddies  are  assumed  to  satisfy 
some  sort  of  balanced  dynamics,  such  as  geostrophy. 
These  balanced  motions  are  also  called  the  “slow  mani- 
fold”. Balanced  turbulent  flows  exhibit  an  upscale  energy 
cascade  and  are  very  inefficient  in  dissipating  energy.  Inter- 
nal gravity  waves  and  other  more  nearly  isotropic  flows,  on 
the  other  hand,  exhibit  a downscale  cascade  with  effective 
dissipation.  The  structure  of  the  equations  governing  bal- 
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anced  flows  is  such  that  they  can  be  integrated  forward  in 
time  only  if  certain  conditions  are  met.  Analysis  of  several 
classical  instability  problems  shows  that  the  fluid- 
dynamical  transition  occurs  in  the  neighborhood  of  these 
integrability  conditions.  It  is  hypothesized  that  violation  of 
these  conditions  characterizes  in  general  the  flow  regimes 
where  energy  is  transferred  from  balanced  to  unbalanced 
motions. 

The  interplay  of  isopycnal  and  diapycnal 
processes 

When  a scalar  is  stirred  by  high  Reynolds  number  tur- 
bulence in  homogeneous  water,  the  eddy  flux  is  mainly 
determined  by  the  turbulent  stirring  and  is  largely  inde- 
pendent of  the  small  value  of  molecular  diffusivity.  Reduc- 
ing the  latter  would  merely  lead  to  finer  scales  in  the  tracer 
field,  without  significantly  affecting  the  turbulent  trans- 
port. The  situation  is  more  subtle  in  the  ocean  because  of 
the  need  to  quantify  both  isopycnal  and  diapycnal  trans- 
ports. One  issue,  already  raised,  is  the  extent  to  which  en- 
ergy released  by  mesoscale  eddies  is  dissipated  in  the 
ocean  interior,  perhaps  by  breaking  internal  waves  which 
have  drawn  some  of  that  energy  from  the  eddies.  This 
could  partly  tie  diapycnal  mixing  rates  to  the  lateral  stir- 
ring rates. 

Another  issue  concerns  the  ultimate  mechanism  for  the 
destruction  of  fine-scale  lateral  gradients  of  temperature 
and  salinity  created  by  stirring  on  isopycnals.  One  possibil- 
ity is  that  vertical  shear  of  the  stirring  process  leads  to  very 
strongly  tilted  frontal  regions,  so  that  pre-existing  diapy- 
cnal mixing  can  remove  the  strong  gradients.  This  could  be 
described  as  a “passive”  destruction  of  the  strong  lateral 
gradients.  A more  active  process  is  possible:  thermohaline 
fronts  caused  by  lateral  stirring  on  isopycnals  could  be- 
come unstable  to  a double-diffusive  interleaving  process. 
This  would  introduce  further  diapycnal  mixing  at  the  same 
time  as  it  destroyed  lateral  gradients. 

The  relative  importance  of  these  scenarios  has  not  been 
established.  Preliminary  estimates  favor  the  latter  process, 
the  signature  of  which  should  be  internal  double-diffusive 
interleaving  regions  in  the  ocean.  These  are  certainly 
known  to  occur,  but  a systematic  catalog  of  their  locations 
and  characteristics  has  not  been  developed.  If  the  process 
is  important  it  may  provide  significant  diapycnal  mixing  in 
regions  with  strong  lateral  gradients  of  potential  tempera- 
ture, and  compensating  salinity,  on  isopycnals.  Moreover, 
the  diapycnal  mixing  is  upgradient  for  density,  as  is  neces- 
sary in  any  double-diffusive  process. 

A significant  sub-theme  which  emerged  on  this  and 
other  issues  concerns  the  extent  to  which  different  stirring 
and  mixing  processes  are  in  parallel  or  in  series.  If  two 


processes  are  in  parallel,  their  strengths  may  be  compared 
to  decide  which  is  more  important.  If  processes  are  in  se- 
ries, one  might  need  to  know  which  is  the  rate-controlling 
one.  In  the  usual  case  of  the  stirring  of  a tracer  in  a homo- 
geneous fluid,  molecular  mixing  clearly  operates  in  series, 
but  with  the  ultimate  mixing  rate  being  controlled  by  the 
turbulent  stirring.  In  the  ocean,  where  there  might  be  some 
intermediate  process  between  the  stirring  agents  and  the 
final  molecular  diffusion,  the  answer  may  be  more  compli- 
cated. For  example,  the  flux  of  warm  water  across  the  Gulf 
Stream  may  be  controlled  by  the  amount  of  air-sea  interac- 
tion that  occurs  in  warm-core  rings  before  they  are  reab- 
sorbed by  the  Stream;  the  flux  is  not  necessarily  controlled 
just  by  the  ring  formation  rate.  The  general  issue  of  series 
or  parallel  processes  needs  to  be  borne  in  mind  in  any  at- 
tempt to  focus  on  processes  that  are  important,  and  may  be 
aided  by  a triple  decomposition  into  mean,  eddies  and  tur- 
bulence (Figure  4). 


Eddies 


Dissipation 


Figure  4.  A simple  schematic  of  the  path  of  scalar  concentration 
variance  from  mean  to  dissipation  via  eddies  and  turbulence. 
(Courtesy  of  C.  Garrett). 


Processes  in  the  surface  mixed  layer 

Because  of  intense  vertical  mixing  there,  the  surface 
mixed  layer  of  the  ocean  provides  many  examples  of  the 
interplay  of  stirring  and  mixing  processes.  One  such  proc- 
ess relies  on  the  fact  that  the  density  of  sea  water  is  deter- 
mined by  both  temperature  and  salinity. 

Typical  of  the  surface  mixed  layer  are  compensated 
fronts  that  are  warm  and  salty  on  one  side  and  cold  and 
fresh  on  the  other  side  such  that  the  density  contrast  is 
small.  Figure  5 shows  an  example  along  a longitude  line  in 
the  subtropical  North  Pacific.  The  existence  of  the  fronts 
has  been  rationalized  as  follows  (Young,  1994):  Nonho- 
mogeneous  atmospheric  forcing  and  entrainment  of  ther- 
mocline  waters  generate  horizontal  gradients  in  tempera- 
ture, salinity,  and  density.  Density  fronts  slump  by  the  ac- 
tion of  horizontal  density-driven  shear  currents  which  are 
eventually  arrested  by  vertical  mixing.  This  process  can  be 


FROM  STIRRING  TO  MIXING  IN  A STRA  TIFIED  OCEAN 


Latitude  (degN) 


Figure  5.  Potential  temperature  (red  line),  salinity  (blue  line)  and 
potential  density  (black  line)  from  a horizontal  Sea  Soar  tow  at 
50  m in  the  subtropical  North  Pacific,  at  140°W,  between  28  and 
30°N.  The  depth  is  in  the  middle  of  the  local  mixed  layer.  The 
vertical  axes  for  temperature  and  salinity  are  scaled  by  the  expan- 
sion coefficients  so  that  excursions  of  temperature  and  salinity 
show  the  change  they  imply  for  density.  (Courtesy  of  R.  Ferrari). 

thought  of  as  shear  dispersion.  Compensated  fronts  do  not 
experience  such  shear  dispersion.  The  result  is  that  density 
fronts  diffuse  whereas  compensated  fronts  persist.  The 
ubiquitous  existence  of  compensated  fronts  suggest  that 
horizontal  mixing  in  the  surface  mixed  layer  is  very  much 
larger  for  density  than  for  “spice”,  the  combination  of 
temperature  and  salinity  that  is  “orthogonal”  to  density. 
Indeed,  eddy-resolving  models  that  diffuse  temperature 
and  salinity  and  hence  density  and  spice  with  the  same 
coefficient  are  unable  to  reproduce  the  observed  tendency 
towards  thermohaline  compensation. 

Air-sea  fluxes  and  the  behavior  of  the  surface  mixed 
layer  are  also  involved  in  a thermodynamic  approach  to 
water-mass  formation.  Basically,  water  must  flow  across 
isopycnals  to  balance  surface  buoyancy  loss  or  gain,  and 
the  convergence  of  this  volume  flux  provides  an  estimate 
of  the  subduction  rate.  This  rate  may  be  compared  with  a 
dynamically-based  prediction  from  Ekman  convergence. 
Both  the  thermodynamic  and  dynamic  approaches,  how- 
ever, require  allowance  for  the  role  of  eddies  and  diurnal 
and  seasonal  cycling  of  the  mixed  layer  depth.  The  prob- 
lem is  thus  complicated  by  the  interplay  of  various  stirring 
and  mixing  processes  which  need  to  be  understood. 

The  stirring  and  mixing  of  non-conservative 
tracers 

Remote  sensing  of  the  sea  surface  from  space  has  made 
us  aware  of  beautiful  and  complex  patterns  in  biologically 
active  variables  such  as  chlorophyll  (Figure  6).  One  ques- 
tion is  whether  the  observed  patchiness  is  a consequence  of 
biological  processes  or  merely  of  the  passive  stirring  of 
large-scale  gradients.  Some  clues  come  from  comparing 
the  shapes  of  the  wavenumber  spectra  of  active  and  pas- 
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sive  scalars,  such  as  chlorophyll  and  temperature  respec- 
tively. If  the  former  shows  comparatively  more  variance  at 
small  scales,  this  is  evidence  for  the  importance  of  scalar 
growth  processes. 

Accounting  for  such  patchiness  in  a model  is  particu- 
larly difficult  if  the  non-conservative  processes  are  nonlin- 
ear. In  that  case  it  is  clearly  inadequate  to  proceed  with  a 
model  that  deals  only  with  an  ensemble  average  variable;  a 
Lagrangian  approach  which  follows  the  development  of 
individual  patches  is  necessary.  A very  basic  model  of  this 
kind  is  the  “Brownian  bug  model”.  It  describes  an  ensem- 
ble of  bugs  that  perform  random  walks  and  experience 
random  birth  and  death  processes  along  the  way.  Mathe- 
matically, this  is  a “superprocess.”  A homogeneous  initial 
distribution  spontaneously  develops  patches  and  voids. 


Figure  6.  Chlorophyll  concentration  in  the  Tasmanian  Sea  as 
observed  by  SeaWiFS.  (Courtesy  of  B.L.  Hua). 

Examining  data  for  evidence  of  patchiness  or  intermit- 
tency  can  be  aided  by  examination  of  high  order  structure 
functions.  For  a variable  C(x,  f)  the  qXh  order  structure 
function  is  defined  as  (8(C(1)) q),  where  8 C{1)  is  the  differ- 
ence in  C over  a separation  / and  ()  denotes  an  ensemble 
average.  The  structure  function  exponent  defined  by 
fitting  a power  law  l^'q)  to  the  structure  function,  character- 
izes the  degree  of  intermittency  for  multifractal  processes. 
It  is  a linear  function  of  q for  non-intermittent  processes 
and  becomes  a nonlinear  concave  function  for  intermittent 
processes.  Calculation  of  this  exponent  from  oceanic  ob- 
servations suggests  that,  in  the  scale  range  of  nearly  iso- 
tropic turbulence,  passive  scalars  such  as  temperature  are 
more  intermittent  than  the  advective  velocity  field.  This  is 
a result  also  obtained  in  numerical  simulations  of  3D  ho- 
mogeneous turbulence  (e.g.,  Jimenez,  2000).  Biologically 
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active  tracers  show  intermittency  beyond  that  associated 
with  turbulent  stirring.  The  excess  increases  with  decreas- 
ing turbulence  intensity  (Figure  7).  These  observational 
results  need  to  be  corroborated  and  understood.  Further  to 
this,  the  possibility  of  a relationship  between  biological 
and  physical  variables  can  be  studied  through  examination 
of  the  behavior  of  joint  structure  functions,  though  caution 
is  required  to  avoid  over-interpretation  of  artifacts  intro- 
duced by  particular  sampling  schemes  or  noise. 


q 

Figure  7.  Empirical  values  of  the  structure  function  exponent 
tjiq)  for  temperature  (circles)  and  in  vivo  fluorescence  for  differ- 
ent values  of  the  turbulent  dissipation  rate  e (red  e = 1(T6  m2  s'3; 
orange  e = 5 x 10'6  nr  s'3;  blue  e = 10'5  m2  s'3;  dark  green  e = 5 
xlO'5  m2  s'3;  light  green  e = 1 O'4  nr  s'3).  The  linear  curve  (dotted 
line)  is  expected  for  non-intermittent  turbulence.  The  data  were 
collected  in  the  eastern  English  Channel.  (Courtesy  of  L.  Seu- 
ront). 


Mixing  efficiency 

Much  of  the  diapycnal  mixing  in  the  ocean  is  associated 
with  the  shear  instability  of  internal  waves.  The  subsequent 
overturning  and  stirring  leads  to  finescale  patterns  of  den- 
sity which  then  disappear  through  the  action  of  molecular 
thermal  diffusivity  or  the  molecular  diffusion  of  salt.  The 
turbulent  eddies  are  being  dissipated  by  molecular  viscos- 
ity at  the  same  time  as  they  are  working  against  buoyancy 
forces;  a major  question  concerns  the  relative  fractions  of 
turbulent  kinetic  energy  lost  to  dissipation  at  a rate  e and 
converted  to  mean  potential  energy  at  a rate  KVN2,  where 
Kv  is  the  eddy  diffusivity  and  N the  buoyancy  frequency. 
The  ratio  T = KvN2/e  is  generally  thought  to  be  about  0.2, 
largely  based  on  calibration  against  other  methods  of  esti- 


mating Kv.  There  is  no  reason  to  suppose  that  T is  a univer- 
sal constant,  however,  but  it  has  not  been  established  what 
environmental  parameters  might  determine  its  value. 

Some  evidence  reported  at  the  meeting,  based  largely  on 
a laboratory  experiment  but  with  some  field  support,  sug- 
gests that  T is  less  than  the  value  of  about  0.2  quoted 
above  if  the  turbulence  is  weak  with  a “turbulent  Reynolds 
number”  ReT=  e/v/V2  less  than  about  100.  This  is  not  sur- 
prising: weak  turbulence  dies  away  before  it  has  accom- 
plished much  buoyancy  flux.  What  was  a surprise  was  a 
claim  that  T also  falls  off  again  for  values  of  Rer  greater 
than  about  1,000.  This  dependence  is  hard  to  accept;  it 
may  transpire  that  ReT  in  these  experiments  was  actually  a 
proxy  for  an  input  variable  that  does  not  involve  the  vis- 
cosity v.  One  variable  that  might  lead  to  reduced  F in  the 
ocean  is  the  duration  of  high  shear  compared  with  the 
buoyancy  time  TV"1.  A long-lived  mixing  event  might  pro- 
duce all  its  buoyancy  flux  at  first  and  subsequently  stir 
water  that  is  already  mixed,  dissipating  energy  but  with 
little  further  buoyancy  flux.  More  study  is  required,  with 
high  resolution  3D  numerical  simulations  of  breaking  in- 
ternal wave  fields  becoming  increasingly  feasible  and  in- 
formative. 

The  question  of  mixing  efficiency  is  important  in  any 
assessment  of  the  amount  of  mixing  that  might  be  associ- 
ated with,  for  example,  the  energy  loss  from  the  barotropic 
tide  to  internal  tides.  There  are  also  implications  for  the 
circulation  that  is  diagnosed  from  measured  dissipation 
rates  £ and  the  observed  mean  density  structure.  The  mean 
diapycnal  velocity  w,  may  be  obtained  from  the  mean 
buoyancy  equation  (e.g.,  St.  Laurent  et  al.,  2001) 

w,N2  = 9(Te)/9z 

so  that  the  value  of  w,  can  depend  on  gradients  of  T.  Vor- 
tex stretching  given  by  fdwjdz  then  affects  the  mean  flow 
along  isopycnals. 

Differential  mixing 

The  rapid  decay  of  turbulent  kinetic  energy  at  small 
values  of  the  turbulent  Reynolds  number  Rer  = dvN2  can 
also  lead  to  less  turbulent  transport  of  salt  than  heat,  given 
the  very  small  molecular  diffusivity  of  the  former.  Al- 
though the  effect  is  clearly  established  in  2D  and  3D  nu- 
merical simulations  and  in  laboratory  experiments  that 
mimic  the  breaking  of  internal  waves  in  a stratified  ocean, 
it  is  somewhat  counterintuitive  at  first  sight.  If  the  turbu- 
lent fluxes  were  downgradient  at  all  scales  then  the  more 
diffusive  tracer  would  lose  its  higher  wavenumber  fluctua- 
tions and  would  experience  less  turbulent  transport.  It 
seems,  however,  that  the  fluxes  are  upgradient  at  high 
wavenumbers,  as  a simple  consequence  of  the  secondary 
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instability  and  collapse  of  an  overturning  billow.  The  more 
diffusive  tracer  thus  loses  some  upgradient  flux  and  there- 
fore experiences  more  downgradient  transport. 

The  importance  of  differential  mixing  for  the  ocean  is 
not  yet  known.  It  may  turn  out  that  in  regions  where  mix- 
ing matters  it  is  at  high  enough  ReT  that  the  eddy  diffusivi- 
ties  of  heat  and  salt  are  the  same.  This  is  particularly  likely 
to  be  true  given  the  intermittent  nature  of  mixing  in  the 
ocean  compared  with  the  rather  steady  stirring  in  the  nu- 
merical experiments;  ReT  based  on  e during  a mixing  event 
in  the  ocean  may  be  much  larger  than  that  based  on  the 
average  s. 

Different  rates  of  mixing  for  heat  and  salt  are,  however, 
firmly  established  for  regions  of  the  ocean  where  double- 
diffusive  processes  are  significant.  These  regions  include 
places,  such  as  the  C-SALT  region  near  Barbados 
{Schmitt,  1994),  where  the  vertical  profiles  of  temperature 
and  salinity  are  nearly  compensating  in  density,  and  places 
such  as  the  Arctic  Ocean  where  internal  wave  activity  is 
low.  Even  in  these  regions  there  are  significant  puzzles, 
partly  because  laboratory  results  seem  not  to  be  readily 
applicable,  and  partly  because  we  lack  adequate  theoretical 
guidance.  One  obvious  problem  is  to  explain  the  step  size 
in  thermohaline  staircases.  In  the  diffusive  regime  (with 
warm  salty  water  below  cooler  fresher  water),  some  useful 
predictions  come  from  a local  water  column  model, 
whereas  in  the  salt  fingering  regime  (with  the  opposite 
stratification)  it  seems  that  the  step  size  might  be  set  by 
lateral  intrusive  processes.  The  issue  is  certainly  not  fully 
resolved,  and  the  interfaces  in  the  fingering  regime  also 
seem  more  complex  in  the  ocean  than  in  the  laboratory. 

Importantly,  many  regions  of  the  ocean  which  might 
support  double-diffusive  processes  seem  to  have  these 
disrupted  by  things  like  the  shear  of  internal  waves,  with 
the  breaking  of  these  waves  then  dominating  diapycnal 
transports.  Even  these  regions  can  show  some  evidence  for 
double  diffusive  transports,  however,  so  that  the  phenom- 
ena warrant  continued  attention. 

Summary 

There  are  typically  nearly  1020  spatial  degrees  of  free- 
dom between  the  smallest  scales  being  resolved  by  nu- 
merical models  of  ocean  circulation  and  the  smallest  scales 
of  variability.  The  meeting  was  therefore  being  ambitious 
in  attempting  to  survey  a vast  range  of  scales  and  phenom- 
ena. Some  general,  if  rather  platitudinous,  principles  do 
seem  clear:  (1)  Studying  the  processes  at  work  in  the  ocean 
is  worthwhile  science.  (2)  There  is  a need  not  just  to  un- 
derstand these  processes  but  also  to  parameterize  them  in  a 
form  that  can  be  implemented  in  models  for  the  large- 
scale  behavior  of  the  ocean.  (3)  These  models  should  be 


used  in  sensitivity  tests  to  determine  which  parameteriza- 
tions  need  improvement. 

It  is  not  really  clear  what  issues  have  the  highest  prior- 
ity. After  a sometimes  tumultuous  discussion  the  partici- 
pants arrived  at  the  following  ranking  by  vote: 

1 . What  are  the  constraints  on  eddy  parameterization 
and  how  should  they  be  implemented? 

2.  What  nonlocal  effects  matter  and  how  can  they  be 
parameterized? 

3.  How  uniform  are  stirring  and  mixing  processes?  Are 

there  any  “hotspots,”  particularly  in  association 
with  topographic  features? 

4.  What  are  the  effects  of  diapycnal  boundary 

processes  on  eddy  parameterization? 

5.  What  is  the  parameter  dependence  of  the  mixing 

efficiency  or  flux  Richardson  number? 

6.  How  does  energy  leak  from  the  slow  manifold? 

7.  What  aspects  of  stirring  and  mixing  affect  biological 

processes?  What  do  patterns  of  biological  tracers 
tell  us  about  fluid  dynamics? 

8.  Under  what  circumstances  do  we  need  different 

diffusivities  for  different  tracers? 

9.  What  are  the  effects  of  adiabatic  stirring  on 

diapycnal  mixing? 

10.  Do  subgridscale  processes  provide  stochastic 
forcing? 

The  last  question  addresses  a very  fundamental  aspect 
of  the  Reynolds  decomposition.  A second  fundamental 
aspect  is  that  the  simple  Reynolds  decomposition  ideally 
requires  a spectral  gap,  in  time  or  space,  separating  the 
unresolved  processes  to  be  parameterized  from  the  slowly 
changing  “mean.”  Both  these  issues  were  discussed  at  the 
meeting,  but  without  any  profound  new  insights. 

Stirring  and  mixing  in  a stratified  ocean  is  the  physics 
that  needs  to  be  parameterized  in  ocean  circulation  models. 
Challenging  open  problems  remain  at  all  levels,  from  very 
fundamental  to  highly  applied  aspects. 
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